Proceedings  of  the  Combustion  Institute,  Volume  28,  2000/pp.  1859-1866 


DYNAMIC  RESPONSE  OE  STRAINED  PREMIXED  ELAMES  TO 
EQUIVALENCE  RATIO  GRADIENTS 


Y.  M.  MARZOUK,^  A.  F.  GHONIEM^  and  H.  N.  NAJM^ 
^Department  of  Meehanieal  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139,  USA 
^Combustion  Research  Facility 
Sandia  National  Laboratories 
Livermore,  CA  94551,  USA 


Premixed  flames  encounter  gradients  of  mixture  equivalence  ratio  in  stratified  charge  engines,  lean 
premixed  gas-turbine  engines,  and  a  variety  of  other  applications.  In  cases  for  which  the  scales — spatial  or 
temporal — of  fuel  concentration  gradients  in  the  reactants  are  comparable  to  flame  scales,  changes  in 
burning  rate,  flammability  limits,  and  flame  structure  have  been  observed.  This  paper  uses  an  unsteady 
strained  flame  in  the  stagnation  point  configuration  to  examine  the  effect  of  temporal  gradients  on  com¬ 
bustion  in  a  premixed  methane/air  mixture.  An  inexact  Newton  backtracking  method,  coupled  with  a 
preconditioned  Krylov  subspace  iterative  solver,  was  used  to  improve  the  efficiency  of  the  numerical 
solution  and  expand  its  domain  of  convergence  in  the  presence  of  detailed  chemistry. 

Results  indicate  that  equivalence  ratio  variations  with  timescales  lower  than  10  ms  have  significant  effects 
on  the  burning  process,  including  reaction  zone  broadening,  burning  rate  enhancement,  and  extension  of 
the  flammability  limit  toward  leaner  mixtures.  While  the  temperature  of  a  flame  processing  a  stoichio- 
metric-to-lean  equivalence  ratio  gradient  decreased  slightly  within  the  front  side  of  the  reaction  zone, 
radical  concentrations  remained  elevated  over  the  entire  flame  structure.  These  characteristics  are  linked 
to  a  feature  reminiscent  of  “back-supported”  flames — flames  in  which  a  stream  of  products  resulting  from 
burning  at  higher  equivalence  ratio  is  continuously  supplied  to  lower  equivalence  ratio  reactants.  The 
relevant  feature  is  the  establishment  of  a  positive  temperature  gradient  on  the  products  side  of  the  flame 
which  maintains  the  temperature  high  enough  and  the  radical  concentration  sufficient  to  sustain  combus¬ 
tion  there.  Unsteadiness  in  equivalence  ratio  produces  similar  gradients  within  the  flame  structure,  thus 
compensating  for  the  change  in  temperature  at  the  leading  edge  of  the  reaction  zone  and  accounting  for 
an  observed  “flame  inertia.”  For  sufficiently  large  equivalence  ratio  gradients,  a  flame  starting  in  a  stoi¬ 
chiometric  mixture  can  burn  through  a  very  lean  one  by  taking  advantage  of  this  mechanism. 


Introduction 

Premixed  flames  encounter  equivalence  ratio  gra¬ 
dients  in  many  applications.  Important  examples  in¬ 
clude  stratified  charge,  spark-ignition  engines  [1,2] 
and  lean  premixed  gas-turbine  combustors  [3,4],  In 
the  former,  a  small  volume  of  the  total  mixture  is 
made  rich  close  to  the  spark  plug,  to  guarantee  ig¬ 
nition  and  help  sustain  flame  propagation  in  the 
leaner  bulk  of  the  mixture.  In  the  latter,  pressure 
oscillations  arising  from  coupling  between  the  acous¬ 
tic  field  and  the  combustion  zone  can  lead  to  oscil¬ 
lations  in  air  mass  flow  rate  near  the  fuel  injectors, 
producing  oscillations  in  equivalence  ratio  at  the 
flame.  Moreover,  pulsed  fuel  injection  may  be  used 
to  impose  similar  oscillations  in  mixture  composition 
to  suppress  combustion  instabilities  resulting  from 
alternate  mechanisms.  In  yet  other  applications  for 
which  burning  occurs  through  nominally  premrxed 
reactants,  a  flame  can  still  encounter  pockets  of 


richer  or  leaner  mixture  produced  by  inhomoge¬ 
neous  mixing,  and  thus  propagate  through  gradients 
in  fuel  concentration. 

In  all  the  cases  described  above,  when  the  tur¬ 
bulence  intensity  is  moderate,  the  combustion  pro¬ 
cess  can  be  modeled  as  a  wrinkled  laminar  flame  in 
which  burning  occurs  across  a  strained  laminar  flame 
surface.  This  paper  examines  the  impact  of  equiva¬ 
lence  ratio  variations  on  the  structure  and  burning 
rate  of  such  a  strained  flame. 

Order-of-magnitude  arguments  can  be  used  to 
quantify  single,  dominant  gradients  in  the  equiva¬ 
lence  ratio  for  many  of  the  applications  men¬ 
tioned.  In  the  case  of  combustion  instability,  oscil¬ 
lation  frequencies  are  typically  0(500-1000  Hz),  and 
hence  the  timescale  of  cj)  variation  is  approximately 
1-2  ms.  In  a  spark-ignition  engine,  the  total  com¬ 
bustion  time  is  0(3  ms),  a  small  fraction  of  which  is 
taken  for  the  transition  between  the  stoichiometric 
mixture  near  the  spark  plug  to  the  leaner  mixture  in 
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LAMINAR  PREMIXED  ELAM ES— Unsteady  Effects 


the  bulk  of  the  cylinder  volume.  Under  these  cir¬ 
cumstances,  the  following  questions  arise:  When  the 
timescales  of  (f)  variation  are  of  the  order  of  magni¬ 
tude  of  the  nominal  flame  timescales,  will  the  flames 
respond  in  a  quasi-steady  fashion;  that  is,  will  they 
adjust  almost  instantaneously  to  the  imposed  strain, 
or  will  some  “flame  inertia”  affect  propagation  into 
the  leaner  mixture?  Will  adjustment  to  different  val¬ 
ues  of  4)  introduce  changes  to  the  flame  structure 
that  have  a  significant  impact  on  the  burning  rate 
and  species  distributions?  And  can  the  presence  of 
these  gradients  explain  the  experimentally  observed 
extension  of  flammability  limits  into  leaner  mixtures 
[1,5]? 

In  this  paper,  we  answer  these  questions  using  an 
unsteady,  strained  flame  code  recently  developed 
[6,7],  Because  the  code  considers  a  flame  stabilized 
in  a  stagnation-point  flow,  imposed  variations  in  ({> 
are  modeled  as  temporal  gradients;  these  can  be  ren¬ 
dered  equivalent  to  spatial  gradients  by  a  Lagrangian 
transformation. 


Formulation  and  Numerical  Solution 


The  impact  of  variable  4)  on  burning  under  finite 
strain  was  modeled  with  an  axisymmetric  stagnation- 
point  flame  configuration,  in  which  a  reactants 
stream  impinged  on  a  stream  of  burned  combustion 
products.  This  formulation  has  been  thoroughly  dis¬ 
cussed  elsewhere  [8-11],  so  only  an  outline  is  pre¬ 
sented  here.  One-dimensional  governing  equations 
were  obtained  by  applying  a  boundary  layer  approx¬ 
imation  across  the  flame  and  solving  along  the  stag¬ 
nation  streamline  x  =  0;  here,  y  was  the  coordinate 
normal  to  the  flame  surface.  Letting  the  imposed 
strain  rate  E  define  the  outer  flow  on  the  reactants 
side,  u„  =  £%,«„=  —  Ey,  and  introducing  the  no¬ 
tation  17  =  m/Mco,  V  =  pv,  we  obtained  the  following 
equations  for  species,  energy,  momentum,  and  mass 
conservation,  respectively: 
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Here  Yj.  is  the  mass  fraction  of  species  k,  while  Dj. 
is  the  corresponding  mixture-averaged  diffusion  co¬ 
efficient  and  Wj.  and  icj,  are  the  molar  weight  and 
molar  production  rate,  respectively.  In  the  remaining 
equations,  Cp  is  the  mixture  specific  heat,  1  is  the 
thermal  conductivity,  /ij.  is  the  molar  enthalpy,  p^  is 
the  density  of  the  reactants  mixture,  and  p  is  the 
dynamic  viscosity  of  the  mixture.  The  enthalpy  flux 
term  was  neglected  in  the  energy  equation  while 
thermal  diffusion  velocity  was  neglected  in  each  spe¬ 
cies  transport  equation.  Smooke  [12]  and  others 
have  shown  that  both  effects  are  unimportant  to  the 
laminar  strained  flame. 

In  the  ensuing  analysis,  strain  rate  was  held  con¬ 
stant.  Unsteadiness  entered  in  the  form  of  the  spe¬ 
cies  and  temperature  boundary  conditions  at  ±“; 
these  conditions  define  the  composition  and  tem¬ 
perature  of  the  two  incoming  streams  of  the  stag¬ 
nation-point  flow. 

y  =  ±00;  =  Y^_i„(f),  T  =  TUt)  (5) 

The  continuity  equation  required  only  one  boundary 
condition,  which  specified  zero  velocity  at  the  stag¬ 
nation  point,  V{y  =  0)  =  0.  The  momentum  con¬ 
servation  equation  required  two  boundary  condi¬ 
tions:  At  an  unburned  stream,  u  =  m„,  so  the 
boundary  condition  is  by  definition  1/  =  1.  Setting 
the  spatial  gradients  in  equation  3  to  zero  gives  the 
boundary  condition  on  the  burned  stream;  for  steady 
strain,  it  becomes  Ur,  =  Jpjph,  where  /?b  is  the 
density  of  the  burned  mixture. 

To  accurately  characterize  flame  structure  and  re¬ 
sponse,  detailed  transport  and  chemistry  were  used 
to  evaluate  transport  coefficients  and  reaction  terms 
in  the  model.  Following  the  approach  validated  in 
Ref.  [12],  a  fixed  Lewis  number  Lej.  was  chosen  for 
each  species,  constant  across  the  flame.  The  Prandtl 
number  was  chosen  constant  and  equal  to  0.75, 
while  a  curve  fit  was  used  to  calculate  the  ratio  i/Cp. 
We  used  a  C(i)  kinetic  model  for  methane/air  com¬ 
bustion  given  in  Ref.  [13],  consisting  of  46  reactions 
among  16  species. 

Numerical  solution  of  the  governing  equations  was 
obtained  via  a  fully  implicit  finite  difference  method, 
as  necessitated  by  the  stiffness  of  detailed  kinetics. 
A  first-order  backward  Euler  formulation  was  used; 
a  first-order  upwind  discretization  was  applied  to 
convective  terms,  while  diffusion  terms  were  discre¬ 
tized  to  second-order  accuracy.  These  spatial  dis¬ 
cretizations  were  performed  on  a  non-uniform  adap¬ 
tive  grid,  permitting  a  dynamic  clustering  of 
meshpoints  in  regions  where  spatial  gradients  were 
strong.  We  enforced  a  limit  on  the  maximum  change 
between  neighboring  meshpoints  for  all  dependent 
variables  (T,  17,  V,  Yj.,  iLj.)  and  their  gradients,  thus 
ensuring  that  the  reaction  zone  of  each  species  was 
resolved  throughout  the  integration  time.  Using 
these  criteria,  for  the  present  mechanism,  the  min¬ 
imum  grid  scale  was  0(10  pm). 


FLAME  RESPONSE  TO  EQUIVALENCE  RATIO  GRADIENTS 
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Fig.  1.  Heat  release  rates  for  static  and  dynamic  flames, 
as  a  function  of  instantaneous  reactants-side  equivalence 
ratio,  (j),.  Arrows  denote  cases  shown  in  Figs.  4,  6,  and  8. 


At  each  time  step,  discretization  reduced  the  gov¬ 
erning  equations  to  a  set  of  nonlinear  algebraic  equa¬ 
tions.  We  implemented  a  novel  set  of  numerical 
methods  to  solve  these  equations  efficiently  and  ro¬ 
bustly  [6].  An  inexact  Newton  iteration  [14]  con¬ 
verged  to  the  solution  of  the  nonlinear  equations  but 
avoided  precisely  solving  the  Newton  condition  far 
away  from  a  solution,  when  the  linear  model  of  New¬ 
tons  method  may  be  poor.  This  inexact  Newton 
method  was  coupled  with  a  backtracking  globaliza¬ 
tion  to  improve  its  domain  of  convergence  [15].  So¬ 
lution  of  the  linear  system  at  each  Newton  iteration 
proceeded  via  a  Kiylov  subspace  method,  Bi- 
CGSTAB  (BiConjugate  Gradients  Stabilized) 
[16,17].  Because  detailed  chemistry  renders  this  lin¬ 
ear  system  ill-conditioned,  BiCGSTAB  must  be  ac¬ 
celerated  with  an  incomplete  LU  factorization  pre¬ 
conditioner  [17]. 


Results 

The  unsteady  strained  flame  model  described 
above  was  used  to  examine  the  impact  of  variation 
on  burning  rate,  flame  structure,  and  extinction. 
Since  the  temperature  and  species  mass  fraction 
boundary  conditions  of  the  formulation  were  arbi¬ 
trarily  adjustable,  numerous  schemes  for  imposing 
mixture  variation  on  the  flame  were  possible.  We 
used  to  denote  the  equivalence  ratio  of  the  re¬ 
actants,  injected  at  y  =  —  with  the  temperature 
of  the  reactants  mixture  always  300  K.  The  compo¬ 
sition  and  temperature  of  the  products  stream,  how¬ 
ever,  corresponded  to  the  chemical  equilibrium  state 
of  a  mixture  with  equivalence  ratio  under  adia¬ 
batic,  constant-pressure  conditions. 


Static  Flames  with  and  without  Back  Support 

We  began  with  a  comparison  of  two  steady-state 
flames,  one  with  =  (^  and  the  other  with 
constant  at  1.0.  The  first  flame  was  labeled  the  NaF 
(for  natural  flame)  and  the  second  the  SuF  (for  sup¬ 
ported  flame).  The  objective  of  this  comparison  was 
to  determine,  in  the  absence  of  time-dependent  ef¬ 
fects,  the  impact  that  products  of  a  mixture 

may  have  on  the  burning  of  a  leaner  mixture.  Spe¬ 
cifically,  we  wished  to  examine  the  mechanism  by 
which  the  high-(^  products  mixture  could  accelerate 
the  burning  of  a  lean  mixture  with  (j>^  above  the  flam¬ 
mability  limit  or  sustain  burning  in  a  lean  mixture 
with  <j)j.  below  the  flammability  limit. 

Figure  1  shows  the  dependence  of  the  integrated 
heat  release  rate  Q  on  <j>^  for  a  number  of  strained 
flames.  All  flames  discussed  in  this  paper  were  sub¬ 
ject  to  a  constant  strain  rate  of  300  s“^.  Here  we 
focus  on  the  static  flames  of  Fig.  1,  marked  with 
triangles.  The  static  simulations  give  rise  to  the  fol¬ 
lowing  observations: 

•  For  all  values  of  4>^,  the  SuF  exhibits  a  higher  heat 
release  rate  than  the  NaF,  with  the  difference  in¬ 
creasing  dramatically  as  the  mixtures  becomes 
leaner. 

•  While  the  NaF  is  extinguished  by  4)^  =  0.5,  the 
SuF  continues  to  burn  down  to  the  lowest  equiv¬ 
alence  ratio  examined. 

The  experimental  results  of  Ref.  [5]  support  the  sec¬ 
ond  observation:  no  matter  how  low  the  equivalence 
ratio  of  propane/ air  reactants,  a  shear  layer  between 
a  reactants  stream  and  a  products  stream  maintained 
at  1770  K  showed  a  finite  rate  of  heat  release.  Ra  [1] 
found  that  it  was  not  possible  to  initiate  a  flame  in  a 
stationary  propane/air  mixture  with  <  0.55  using 
a  laser-ignition  source  in  experiments,  or  with  < 
0.4  using  a  heat  source  in  numerical  simulations. 
However,  igniting  a  small  bubble  of  stoichiometric 
mixture  did  sustain  burning  for  some  time  in  a  mix¬ 
ture  with  (f)  slightly  below  0.5  in  the  experiments, 
and  for  even  lower  values  of  (j)  in  numerical  simu¬ 
lations. 

To  examine  the  reasons  for  enhanced  or  sustained 
burning  in  the  SuF,  Fig.  2  shows  profiles  of  tem¬ 
perature,  OH  mass  fraction  and  diffusion  flux,  and 
the  source  term  in  the  energy  equation,  w^  =  Sj. 
wfli]^,  for  the  NaF  and  the  SuF  at  (/),  =  0.6.  Several 
interesting  observations  can  be  made  here: 

•  While  the  temperature  at  the  peak  burn  rate  is 
almost  the  same  for  both  flames,  w^  is  much 
higher  in  the  SuF  than  in  the  NaF. 

•  The  tail  of  liy  on  the  burned  side  of  the  SuF  is 
longer  than  that  of  the  NaF;  that  is,  burning  is 
sustained  further  on  the  products  side  of  SuF,  con¬ 
tributing  to  the  enhancement  in  Q. 

•  As  the  burning  rate  of  a  stagnation-point  flame  de¬ 
creases,  the  reaction  zone  tends  to  stabilize  itself 
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Fig.  2.  Supported  flame  and  natural  flame  structures, 
static  at  d  shosving  T,  tix,  Yoh>  and  net  diffusion 

of  OH. 


Fig.  3.  Q  versus  time  for  natural  flames  subject  to  vari¬ 
ous  d<l>/dt.  Lines  are  dynamic  behavior,  while  the  symbols 
are  static  Q  corresponding  to  the  instantaneously  imposed 
</).  Arrows  denote  cases  shown  in  Figs.  4  and  6. 

closer  to  the  stagnation  plane.  Indeed,  in  the  illus¬ 
trated  case  of  4>^  =  0.6,  the  reaction  zone  of  the 
NaF  straddles  the  stagnation  point  while  the  SuF 
is  stabilized  well  into  the  reactants  stream. 

Since  the  reaction  zone  temperatures  of  the  two 
flames  were  nearly  equal,  a  strong  contribution  to 
the  difference  in  reaction  rates  must  be  the  higher 
concentration  of  radicals  in  the  reaction  zone  of  the 
SuF — for  OH,  nearly  twice  that  of  the  NaF.  Both 
temperature  and  radical  concentrations  imposed  at 
the  products-side  boundary  condition  were  higher 
in  the  SuF  than  in  the  NaF,  thus  maintaining  higher 
values  throughout  the  products  side,  with  radical 
concentrations  increasing  further  to  the  right  of  the 
reaction  zone.  The  result  was  a  stronger  net  diffusive 


flux  of  OH  into  the  reaction  zone  of  the  SuF,  as  seen 
in  Fig.  2.  Similar  results  were  observed  for  the  rad¬ 
icals  O  and  H. 

As  falls  below  0.6,  the  SuF  continues  to  burn 
while  extinction  is  observed  in  the  N aF.  This  precip¬ 
itous  decline  in  reaction  rates  coincides  with  the  re¬ 
action  zone  of  the  NaF  crossing  the  stagnation  plane 
into  regions  of  high  products  concentration.  While 
the  temperature  in  the  reaction  zone  remains  high, 
fuel,  oxidizer,  and  radical  concentrations  drop  due 
to  the  presence  of  products,  leading  to  a  falloff  of  all 
reaction  rates  [18].  Simultaneously,  convection  re¬ 
verses  sign  and  counteracts  diffusion,  driving  reac¬ 
tants  away  from  the  reaction  zone  and  toward  the 
stagnation  plane. 

Static  versus  Dynamic  Natural  Flames 

To  examine  the  impact  of  (f)  variation  on  heat  re¬ 
lease  rate  and  flame  structure,  we  turned  to  cases  in 
which  varies  linearly  in  time.  For  these  cases, 
immediately  ahead  of  the  flame  was  reduced  from 
1.0  to  0.3  over  four  distinct  time  spans:  1.75  ms,  3.5 
ms,  17.5  ms,  and  35  ms.  Here,  “ahead  of  the  flame” 
is  defined  as  the  leading  edge  of  the  temperature 
profile,  the  point  at  which  the  temperature  gradient 
rises  to  a  small  finite  threshold. 

Definition  of  the  mixture  boundary  conditions  on 
the  burned  side  of  the  flame  is  more  subtle,  for  as  a 
strained  flame  propagates  through  an  equivalence 
ratio  gradient,  the  state  of  the  products  left  in  the 
wake  of  the  flame  may  change  continuously.  To 
model  this  case  using  the  stagnation-point  flame 
code,  we  note  that  the  mixture  takes  essentially  a 
flame  timescale  to  burn  through  the  flame  structure, 
that  is,  to  reach  the  products  side.  Therefore,  if  the 
conditions  on  the  burned  side  are  determined  by  the 
combustion  of  the  reactants  as  they  move  through 
the  flame,  and  if  the  reactants  state  changes  contin¬ 
uously,  conditions  on  the  burned  side  must  be  de¬ 
layed  from  conditions  on  the  reactants  side  by  a 
flame  timescale.  Thus,  we  assumed  that  </>(,  =  <^^(1 
—  Tf)  and  updated  both  continuously,  taking  a  nom¬ 
inal  flame  time  scale  Tf  to  be  1  ms.  A  flame  using  this 
scheme  for  defining  boundary  conditions  was  la¬ 
beled  a  DNaF  {Dynamic  Natural  Flame). 

To  assess  the  impact  of  variation,  we  compared 
heat  release  rates  for  these  dynamic  flames  with 
those  obtained  for  static  flames,  shown  as  functions 
of  time  in  Fig.  3.  Static  heat  release  rates  correspond 
to  a  NaF  burning  a  steady  reactants  stream  whose 
4>^  is  the  same  as  the  instantaneous  4>^  ahead  of  the 
dynamic  flame.  The  same  comparison,  plotted 
against  the  instantaneous  value  of  <^)j,  is  shown  in  Fig 
1.  At  the  early  stages,  the  burning  rate  lags  behind 
the  value  corresponding  to  its  current  by  about  1 
ms.  This  is  the  flame  timescale,  the  time  it  takes  the 
reactants  ahead  of  the  flame  to  diffuse  into  the  re¬ 
action  zone  and  burn,  suggesting  that  the  observed 
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Fig.  4.  Comparison  of  dynamic  and  static  natural  flame 
structure  at  Q  =  2.96  X  10®  erg/(cm^  s). 


Fig.  5.  T  at  maxlti;-]-)  versus  Q,  dynamic  and  static  flames. 


time  lag  is  due  to  the  finite  transport  rate  within  the 
flame  structure.  (Note  the  distinction  between  the 
delay  time,  used  to  define  the  boundary  condition 
on  the  DNaF,  and  the  time  lag,  observed  in  flame 
response  to  imposed  changes  in  (/>.)  The  observed 
time  lag  increases  as  decreases  since  the  flame 
becomes  thicker  as  the  mixture  becomes  leaner.  Re¬ 
sults  in  Figs.  1  and  3  indicate  that 

•  the  value  of  Q  is  always  higher  for  the  dynamic 
flame  than  for  the  static  flame,  since  (l>r  is  decreas¬ 
ing. 

•  while  the  static  flame  tends  to  extinguish  around 
<l>^  =  0.5,  the  dynamic  flame  continues  to  burn 
through  lower  values  of  4>r- 

•  the  higher  d(j)Jdt,  the  lower  the  equivalence  ratio 
which  can  support  burning;  that  is,  there  is  a 
“flame  inertia”  whose  effect  becomes  more  pro¬ 
nounced  with  faster  gradients  in  (j>. 

While  the  transport  time  lag  can  explain  the  higher 
burning  rate  of  the  DNaF,  the  interaction  of  finite 


transport  rates  and  d4>ddt  must  have  an  effect  on 
flame  structure.  In  Fig.  4,  we  compare  the  tempera¬ 
ture,  OH  concentration,  and  iv-y  profiles  of  a  dy¬ 
namic  flame  with  those  of  a  static  flame  that  pro¬ 
duces  the  same  value  of  Q.  Both  flames  are  marked 
with  arrows  in  Figs.  1  and  3.  The  static  flame  cor¬ 
responds  to  =  0.7,  and  the  dynamic  flame  cor¬ 
responds  to  the  1.75  ms  change  in  probed  2.075 
ms  from  the  start  of  the  calculations.  The  reason  for 
choosing  the  two  flames  at  the  same  value  of  Q  is  to 
remove  the  impact  of  time  lag  and  different  values 
of  (/)j,  and  to  focus  instead  on  the  effect  of  the  (p 
gradients.  Several  observations  can  be  made: 

•  The  dynamic  flame  stands  further  away  from  the 
stagnation  point;  it  has  a  higher  displacement  ve¬ 
locity  than  the  static  flame. 

•  The  static  flame  has  a  higher  peak  burn  rate  Wr, 
while  the  dynamic  flame  shows  a  wider  reaction 
zone  with  a  longer  tail  on  the  products  side. 

•  The  temperature  of  the  static  flame  is  higher  than 
that  of  the  dynamic  flame  throughout  the  reaction 
zone,  in  particular,  at  the  point  where  Wy  takes  its 
maximum. 

•  OH  concentrations  are  almost  the  same  on  the  re¬ 
actants  side  of  both  flames,  but  higher  on  the  prod¬ 
ucts  side  of  the  dynamic  flame. 

These  distinctions  persist  throughout  the  calcula¬ 
tion.  The  temperature  at  max(!i;i-)  is  shown  for  both 
static  and  dynamic  flames  in  Fig.  5;  at  any  given  Q, 
a  larger  dp/dt  produces  lower  reaction  zone  tem¬ 
peratures.  Clearly,  the  unsteady  change  in  influ¬ 
ences  the  flame  structure  significantly,  and  one  can¬ 
not  explain  Q  enhancement  in  the  dynamic  case 
using  the  transport  lag  argument  alone. 

For  the  same  Q,  lower  temperature  and  reaction 
rates  within  the  dynamic  flame  reaction  zone  are  bal¬ 
anced  by  elevated  temperature  and  OH  concentra¬ 
tion  on  the  products  side;  these  are  the  most  impor¬ 
tant  features  of  Fig.  4.  Similar  to  the  elevated 
products-side  temperature  and  radical  concentra¬ 
tions  introduced  externally  in  the  static  SuF,  these 
gradients  in  the  wake  of  the  reaction  zone  maintain 
higher  wy  in  the  tail  of  the  DNaF.  The  effect  is  seen 
even  more  clearly  in  Fig.  6,  which  shows  the  same 
two  cases  of  static  and  dynamic  flame  structure  as 
Fig.  4.  The  CO  concentration  remains  significantly 
higher  on  the  products  side  of  the  reaction  zone  for 
the  dynamic  flame  than  for  the  static,  before  drop¬ 
ping  to  its  equilibrium  value.  This  CO  tail  corre¬ 
sponds  to  ongoing  CO  oxidation  on  the  products  side 
of  the  flame  and  thus  the  persistence  of  heat  release. 
In  cases  of  lower  d(l>Jdt,  these  products-side  gradi¬ 
ents  in  radical  concentration,  CO  concentration,  and 
temperature  are  not  observed. 

Changes  in  flame  structure  resulting  from  un¬ 
steadiness  in  (j>  were  expected,  since  the  dynamic 
flame  processes  a  mixture  with  a  fuel  concentration 
gradient.  Indeed,  proceeding  from  the  leading  edge 
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Fig.  6.  Static  and  dynamic  CO  profiles,  for  natural 
flames  at  (i  =  2.96  X  10®  erg/(cm^  s)  and  supported 
flames  at  (i  =  5.24  X  10®  erg/(cm^  s). 


Fig.  7.  Q  versus  time  for  supported  flames  subject  to 
various  d4>,/dt.  Lines  are  dynamic  behavior,  while  the  sym¬ 
bols  are  static  Q  corresponding  to  the  instantaneously  im¬ 
posed  (j>r.  Arrows  denote  cases  shown  in  Figs.  6  and  8. 

of  the  flame  through  the  products  side,  we  observed 
the  results  of  burning  at  earlier,  higher  values  of  — 
elevated  temperature,  radical  concentration,  and 
CO  concentrations.  This  “history”  embedded  within 
the  flame  structure  is  a  mechanism  for  the  observed 
“flame  inertia”  of  Fig.  3. 

Static  versus  Dynamic  Supported  Flames 

To  refine  the  results  of  the  previous  section,  we 
examined  the  impact  of  nonzero  d(j)ddt  while  hold¬ 
ing  constant  and  equal  to  1.0,  calling  an  unsteady 
flame  with  this  products-side  boundary  condition  a 
DSuF  (for  dynamic  supported  flame) .  The  DSuF  is 
interesting  for  several  reasons.  First,  it  eliminates  a 
lingering  ambiguity  in  the  structural  comparison  of 


Fig.  8.  Comparison  of  dynamic  and  static  supported 
flame  structure  at  (j)  =  5.24  X  10®  erg/(cm^  s). 


dynamic  and  static  NaFs:  static  and  dynamic  NaFs 
at  the  same  Q  may  not  have  identical  products-side 
boundary  conditions,  due  to  temporal  variation  in 
<^,  while  static  and  dynamic  SuFs  always  have  iden¬ 
tical  products  streams.  The  DSuF  has  physical  rele¬ 
vance  as  well;  the  burning  of  a  lean  pocket  in  an 
inhomogeneously  mixed  reacting  flow,  or  the  burn¬ 
ing  of  a  premixed  flame  on  a  shear  layer,  are  both 
supported  by  a  reservoir  of  hot  combustion  prod¬ 
ucts. 

The  same  linear  changes  in  applied  to  the 
DNaF  were  repeated  for  the  dynamic  supported 
flame.  Fig.  7  shows  Q  versus  time  for  the  four  dy¬ 
namic  cases,  along  with  the  corresponding  static  val¬ 
ues.  The  same  results  are  shown  in  Fig.  1,  plotted 
versus  the  instantaneous  value  of  As  with  the 
DNaF,  the  burning  rate  for  the  dynamic  flame  is 
always  higher  than  that  of  the  corresponding  static 
case.  In  comparison  with  the  DNaF,  however,  we 
observed  the  following: 

•  The  time  lag  was  slightly  smaller  in  the  DSuF  since 
the  flames  were  thinner. 

•  The  DSuF  sustained  burning  at  lower  (f)^,  even 
when  the  (j)  gradients  were  weak,  due  to  the  same 
back  support  discussed  for  the  static  case. 

•  For  fast  changes  in  <t>^,  the  DSuF  and  DNaF 
burned  at  nearly  the  same  rates,  while  for  slower 
changes,  the  burning  rate  of  the  DSuF  was  higher 
than  that  of  the  DNaF,  again  showing  the  effect  of 
back  support. 

To  examine  the  role  of  flame  structure  in  this  dy¬ 
namic  response,  we  again  compared  static  and  dy¬ 
namic  flame  structures  at  equal  Q,  in  Fig.  8.  The 
static  supported  flame  corresponds  to  =  0.8, 
while  the  dynamic  flame  corresponds  to  the  1.75  ms 
change  in  probed  1.55  ms  from  the  start  of  cal¬ 
culations.  As  in  the  comparison  for  natural  flames, 
the  static  flame  sits  closer  to  the  stagnation  point. 
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has  a  higher  peak  burn  rate,  and  has  a  higher  tem¬ 
perature  throughout  the  reaction  zone.  Fig.  5  con¬ 
firms  that  larger  produces  lower  temperatures 
at  the  peak  of  for  any  given  Q.  The  dynamic  flame 
compensates  for  a  lower  peak  in  ivj-  with  a  wider 
reaction  zone — a  longer  tail  in  on  the  products 
side.  This  tail  is  supported  by  higher  temperature 
and  radical  concentrations  on  the  products  side,  and, 
referring  to  Fig.  6,  an  elevated  level  of  CO,  indicat¬ 
ing  ongoing  CO  oxidation.  Because  the  dynamic  and 
static  supported  flames  have  precisely  the  same  in¬ 
coming  products-stream  composition,  differences  in 
flame  structure  to  the  right  of  the  reaction  zone  can 
only  be  due  to  unsteadiness  in  the  reactant  mixture. 
The  equivalence  ratio  gradient  thus  creates  internal 
gradients  in  the  wake  of  the  reaction  zone  regardless 
of  the  products  stream  boundary  condition. 


Conclusions 

Strained  premrxed  flames  propagating  through 
gradients  in  equivalence  ratio  in  which  the  length/ 
timescales  are  comparable  with  the  flame  length/ 
timescale  exhibit  a  substantially  different  response 
than  those  burning  in  uniform  streams.  They  burn 
faster  and  can  burn  into  mixtures  whose  equivalence 
ratios  are  lower  than  those  corresponding  to  flam¬ 
mability  limits.  The  mechanism  responsible  for  these 
changes  is  the  establishment  of  spatial  gradients  in 
temperature  and  the  accompanying  radical  concen¬ 
trations  in  the  wake  of  the  flame — gradients  that  act 
as  a  temporary  “heat  (and  chemical)  reservoir”  for 
new  incoming  reactants.  As  a  result,  for  the  same 
heat  release  rate,  the  reaction  zone  of  a  flame  burn¬ 
ing  in  a  stoichiometric-to-lean  equivalence  ratio  gra¬ 
dient  is  broader  than  that  of  the  nominal  static  flame, 
with  a  lower  temperature  at  the  peak  and  a  higher 
radical  concentration  throughout.  These  effects  are 
associated  with  the  burning  of  a  mixture  that 
changes  its  fuel  content  over  the  thickness  of  the 
flame,  here  falling  to  a  leaner,  more  weakly  burning 
mixture  as  one  proceeds  from  the  products  side  to 
the  reactants  side  of  the  flame.  In  this  sense,  a 
strained  flame  burning  a  sufficiently  steep  equiva¬ 
lence  ratio  gradient  generates  its  own  back  support, 
regardless  of  the  actual  composition  of  its  products 
stream. 

It  is  interesting  to  note — or  at  this  point,  specu¬ 
late — that  a  lower  temperature  within  the  reaction 
zone  resulting  from  an  equivalence  ratio  gradient 
would  lead  to  a  lower  NO^-  formation  rate  for  the 
same  overall  burning  rate.  This  does  not  come  at  the 
price  of  higher  CO,  since  the  temperature  and  the 
radical  concentrations  remain  high  for  CO  oxidation 
on  the  products  side  of  the  flame.  This  may  help 
explain  why  periodic  oscillations  in  equivalence  ratio 
have  been  observed  to  reduce  overall  emissions  [19]. 
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COMMENTS 


Robert  W.  Pitz,  Vanderbilt  University,  USA.  These  weak 
limit  premixed  flames  supported  by  hot  products  are  very 
interesting  because  of  their  relevance  to  direct-injection 
spark-ignition  engines  where  such  flames  are  likely  to  exist. 
An  interesting  question  is  how  much  additional  burning  of 
the  weak  premixed  reactants  occurs  due  to  the  presence  of 
the  hot  products  stream.  In  your  work,  you  calculated  the 
integrated  heat  release  rate  that  is  dominated  by  the  hot 
product  stream.  Did  you  compare  your  integrated  heat  re¬ 
lease  values  for  weak  reactants  vs.  hot  products  as  opposed 
to  hot  products  vs.  air?  The  difference  between  these  two 
would  give  the  additional  heat  release  due  to  burning  of 
the  weak  premixed  mixture  from  the  presence  of  the  hot 
products  stream. 

Author’s  Reply.  We  do  not  expect  significant  heat  release 
to  result  from  air  at  300  K  impinging  on  hot  products.  This 


is  because  the  lower  values  of  the  temperature  may  quench 
CO  oxidation  in  the  hot  products  stream.  Consider  the  CO 
profiles  for  static  flames  in  Fig.  6:  proceeding  from  the 
products  side  to  the  reaction  zone.  CO  concentrations  rise 
monotonically  from  their  equilibrium  values  in  the  burned 
products  mixture,  that  is,  the  CO  concentration  every¬ 
where  stays  at  or  above  its  equilibrium  value  corresponding 
to  the  products  stream  equivalence  ratio.  Comparing  the 
integrated  heat  release  rates  of  supported  and  natural 
flames  thus  provides  a  meaningful  measure  of  how  much 
additional  bnrning  of  the  “weak”  premixed  reactants  occnrs 
due  to  the  presence  of  the  hot  (fif  =  1.0)  prodncts  stream. 
The  higher  heat  release  rate  of  supported  flames  burning 
weak  reactants  is  due  to  the  diffusion  of  heat  and  radicals 
from  the  products  stream. 


